In-situ induced polarization method for determining formation permeability

ABSTRACT

An in-situ method of determining the pore size distribution, capillary pressure curve, and permeability of a formation using induced polarization logging measurements. The induced polarization (IP) voltage decay curve is measured opposite the formation of interest and the voltage decay curve is decomposed into a series of relaxation times related to the pore size distribution. The capillary pressure curve and the permeability are approximated from the pore size distribution. The measurements can also be made in the frequency domain using various frequencies to obtain the variation in phase shift with frequency. The variation in phase shift with frequency can be correlated with previous core measurements to obtain the pore size distribution, capillary pressure curve and formation permeability.

This is a division, of application Ser. No. 591,140, filed Mar. 19, 1984now U.S. Pat. No. 4,644,283, issued Feb. 17, 1987.

BACKGROUND OF THE INVENTION

The present invention relates to an in-situ method for determining thepore size distribution, capillary pressure curve and permeability of aformation surrounding a well. These are all important considerations indeciding whether hydrocarbon-bearing formations are commercial. Aformation may contain a large amount of hydrocarbon but if the pores aretoo small and the permeability is too low it may not be possible toproduce the hydrocarbons commercially. Thus, the determination of thepore size distribution, capillary pressure curve, and permeability ofthe formation is an important consideration in determining whether awell which has penetrated a hydrocarbon-bearing formation should becompleted. The expense of completing a well is, of course, considerablesince the well must first be cased and production tubing and well-headequipment installed.

At present there are only two reliable methods for determining the poresize distribution, capillary pressure curve and permeability of aformation. The method used most often is to take cores from theformation during the drilling of the well and then to analyze the coresin a laboratory. Coring enables high precision laboratory measurementsof capillary pressure and permeability. Coring, however, is veryexpensive and involves some risk. While coring a well, one does not haveas much control over the well as when drilling with conventional rockbits and the possibility of sticking the drill string in the well isgreatly increased during coring operations. Further, during coringoperations, the penetration rates are greatly reduced, thus increasingthe time required to drill the well to a target depth. In addition,measurements on core plugs may not be a good representation of theformation as a whole, since core plugs are small, discrete samples.

The only reliable logging method for determining pore size distribution,capillary pressure curves and permeability is the use of the nuclearmagnetism logging (NML) tool. In this method the T₁ decay curve ofprotons magnetically polarized in the formation is recorded and thenmathematically inverted to obtain a pore size distribution and anapproximate permeability (J. D. Loren and J. D. Robinson, "RelationsBetween Pore Size, Fluid and Matrix Properties, and NML Measurements",Society of Petroleum Engineers Journal, September 1970, pages 268-278).The use of nuclear magnetism logging is commercially practiced andreadily available. While the technique is available, it does haveseveral disadvantages. Among the disadvantages are the small depth ofinvestigation of the nuclear magnetism log (about 6 inches) and therequirement that the borehole fluid be treated with magnetized particlesto eliminate the borehole mud response from the log. Normally, magnetiteparticles in suspension are added to the drilling mud to suppress theborehole response. This involves considerable expense and effort,because it is necessary to circulate the magnetite in the well to assureadequate mixing, and the circulated mud must be checked to verify thatthe nuclear magnetic response has been eliminated.

Another disadvantage of the nuclear magnetism log is that the signal isextremely weak, often requiring several measurements at a fixed locationfor signal averaging. Thus, a good NMR decay curve of the qualityrequired for capillary pressure and permeability determination isusually not obtained while continuously logging. Another disadvantage ofthe NML is that the nuclear magnetism decay times are so short in someformations, for example in tight gas sands, that the signals cannot bemeasured by existing NML tools (J. A. Brown, L. F. Brown, J. A. Jackson,J. V. Milewski, B. J. Travis, "NMR Logging Tool Developments: Laboratorystudies of Tight Gas Sands and Artificial Porous Material", SPE/DOE10813, pages 203-208). This is a severe disadvantage because it isprecisely the marginal, low permeability formations where accuratemeasurement of the permeability is required most. Still anotherdisadvantage of the NML is that the permeability determined from nuclearmagnetism is a three-dimensional average, because the polarized protonson the water molecules diffuse randomly until they relax at the porewalls. In many formations there is an order of magnitude differencebetween vertical and horizontal permeability. The NML results will yielda three-dimensional averaged permeability rather than a separatevertical and horizontal permeability. Yet another difficulty withnuclear magnetism logging is that in oil-bearing formations, both theprotons in water and in oil contribute to the T₁ decay curve. Since thewater phase and the oil phase have different decay rates, an unambiguousdetermination of pore size often cannot be made. In addition, the T₁decay times depend on the character of the pore surface, such ashydrocarbon wetting, bound hydration ions and paramagnetic centers. Asis well known in the art, these can all produce significant effects onthe measured T₁ relaxation times (J. A. Glasel, "NMR Relaxation inHeterogeneous Systems", Nature 227, 704-705 (1970); R. J. S. Brown andI. Fatt, "Measurements of Fractional Wettability of Oilfield Rocks bythe Nuclear Magnetic Relaxation Method", Pet. Trans. AIME 207, 262-264(1956).)

A final problem in relating the NMR relaxation times to formationpermeability is that the pores probed by NMR need not be hydraulicallyconnected. Therefore an impermeable medium containing disconnected vugscould yield the same T₁ decay curve as a permeable rock containingconnected pores.

SUMMARY OF THE INVENTION

The present invention provides a method using in-situ logging methodsfor measurement of the capillary pressure curve and formationpermeability. The method utilizes an induced polarization logging tool,preferably a focused induced polarization logging apparatus and methodsuch as described in copending application Ser. No. 505,623 filed onJune 20, 1983 now U.S. Pat. No. 4,583,046 issued Apr. 15, 1986. In thispatent application there is described an induced polarization loggingtool having means for focusing the electrodes to provide eitherhorizontal or vertical measurements of the induced polarization.Measurements of the logging tool are recorded at the surface in the formof time domain IP decay curves or in the form of phase shift versusfrequency of the applied electric field.

The present invention has several advantages compared to the NML tooldescribed previously. The focused electric field in this inventionpenetrates deep into the formation away from the borehole. In addition,the present method does not require any treatment of the boreholefluids. A major advantage is that the signal to noise in the presentinvention is sufficient that signal averaging is not required and a wellcan be logged continuously. Another advantage of the present inventionis that the IP voltage decay times are long enough even in low marginalreservoirs. Still another advantage is that the ions polarized by thepresent invention move along the direction of the applied electricfield, so that both horizontal and vertical permeability can beseparately determined. Yet another advantage is that, unlike the NML,the IP log observes only the water and not the oil in the formation.This is because oil is substantially nonconducting. A final advantage isthat, unlike the nuclear magnetic response, the IP response is measuredonly along conducting pathways through the formation; therefore,disconnected and dead-end pores, which do not contribute topermeability, are not measured by the IP logging tool.

The method of the present invention requires that the above mentionedlogging tool be operated opposite the formation of interest to measurethe time-domain IP decay curve. The time domain IP decay curve isdecomposed by a computer algorithm into a set of exponential decays withdifferent time constants. From the amplitude and time constant of eachexponential, a pore size histogram is constructed. In particular, therelationship between diffusion constant D of sodium chloride electrolytein aqueous solution (cm² /sec), the displacement length r (cm), and thetime constant t (sec) is used in the expression: ##EQU1##

The value of the diffusion coefficient at the formation temperature isused in this expression. The diffusion constant of sodium chlorideelectrolyte is known to be 1.5×10⁻⁵ cm² /sec at 25° C., and the effectof temperature on the diffusion constant is well known. We havediscovered that the displacement length r is substantially the same asthe pore size in clastic rocks. Thus, the IP decay curve yields thedistribution of pore sizes in the formation.

Once the distribution of pore sizes is known, an approximate value forformation permeability can be computed using various geometric models,which are well known in the art. By proper calibration with core plugsfrom the formation of interest, the response of the induced polarizationlogging tool can be made directly in permeability units.

An alternative method of the present invention is to operate the IPlogging tool of the aforementioned patent in the frequency domain. Thisrequires that the logging tool be operated at several substantiallydifferent frequencies adjacent to the formation of interest. The phaseshift for each frequency is measured and plotted to obtain a curverepresenting the phase shift versus frequency over a substantial rangeof frequencies. This plot can also be obtained from the time domain IPdecay curve described above using linear transform techniques well knownin the art. The phase shift versus frequency curve can be decomposedinto a set of amplitudes and time constants which are used with theabove mentioned relation between diffusion constant D and displacementlength r to obtain a pore size histogram.

A shortened form of the above procedure is to obtain the phase shiftversus frequency plot as described above, and record the frequency(f_(min)) at which the phase shift starts to decrease rapidly towardzero with the frequency dependence substantialy as f¹.0. This frequencyrepresents the largest pore size in the formation, again obtained fromequation (1). Since the largest pores control the formationpermeability, an approximate value of permeability can be obtained fromthe measurement of f_(min) alone.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will be more easily understood from the followingdescription when taken in conjunction with the drawings in which:

FIG. 1 represents an idealized pore in a shaly sand and the chargescarried by clay counterions and by sodium and chloride ions in the freeelectrolyte.

FIG. 2 represents a plot of the induced polarization decay voltagepresent in the formation as the induced polarization logging tool isoperated.

FIG. 3 represents the decomposition of the normalized inducedpolarization decay into a weighted series of four exponential decays.

FIG. 4 represents the cumulative pore volume versus pore size curvecalculated from the induced polarization voltage decay.

FIG. 5 represents a Hg/air and oil/water capillary pressure curvedetermined from the induced polarization voltage decay, compared withthe measured Hg/air capillary pressure curve.

FIG. 6 represents the IP decay curves for the Berea sandstone 100%water-saturated and 50% oil-saturated.

FIG. 7 represents a series of curves of phase shift versus frequency forvarious formations.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

According to the teachings of this invention, there is shown in FIG. 1 aschematic of a pore contained within a shaly sand. The pore is shown tobe spherical in shape with radius A, and has a plurality of pore throatsof radius R less than A which enter the pore. The pore contains a seriesof clay-rich and clay-free zones. The induced polarization logging toolof the aforementioned patent applies a substantially constant electricfield to the formation in the direction shown in FIG. 1. The clay-richzones act as cation-selective membranes and restrict the flow ofnegatively charged cations under the influence of the electric field.Thus, after sufficient time an electrochemical gradient is establishedacross the clay-free zones by the buildup of electrolyte concentrationat the edge of the clay-rich zone. When the electric field from theinduced polarization logging tool is turned off, the concentrationgradient disappears as the electrolyte ions diffuse back to theirequilibrium positions. This, in turn, decreases the induced potentialgradient which constitutes the induced polarization decay curve.

A typical time-domain IP decay curve is illustrated in FIG. 2. Thevoltage Vp is the equilibrium voltage after the electric field has beenapplied for a substantial time, and the voltage V_(s) (t) is the timedomain IP decay voltage as a function of time. The voltage V_(s) (t=0)is the IP decay voltage at the instant the electric field is terminated.

If the porous medium consisted of clay-free zones all of the same size,the induced polarization decay curve would consist of a singleexponential decay of time constant, t: ##EQU2## where D is the diffusionconstant of electrolyte ions in the aqueous solution at formationtemperature, and r is the length of the clay-free zone along thedirection of the applied electric field. For sodium chloride solution at25° C., D=1.5×10⁻⁵ cm² /sec; r is in centimeters, and t is in seconds.However, naturally-occurring sandstones have a wide range of clay-freedistances, so that the observed induced polarization decay curve is acomplex exponential decay consisting of a weighted series ofexponentials with different time constants. The decomposition of the IPdecay curve of FIG. 2 into four exponentials is shown in FIG. 3.

The longest time constant in the IP decay curve will be contributed bythose ions with the longest relaxation distance. For times longer thanthe longest time constant the ions have enough time to return to theirequilibrium positions and the IP decay voltage will have decreased tozero.

We have discovered that the clay-free relaxation distances aresubstantially the same as the pore size distribution. The explanationfor this is that the clay occurs naturally as a pore lining material. Inparticular, referring to FIG. 1, the clay particles in the pore throatare the most effective membranes for the electrolyte contained in thepore. Thus, the scale of relaxation distances is substantially the sameas the pore scale. This discovery has not been previously recognized inthe art and is the basis for determining a pore size histogram from aninduced polarization logging tool.

The preferred method of this invention is to measure an IP decay curveV_(s) (t) opposite the formation using the focused IP logging tool orcopending patent application. Although a focused induced polarizationlogging tool is the preferred embodiment, the method of this inventioncan also be applied to conventional induced polarization logging tools.The IP logging tool applies a substantially constant current into theformation until equilibrium is substantially reached. This time istypically in the range 0.01-100 seconds. The current is then terminatedand the IP decay curve is measured. The IP decay curve V_(s) (t) isnormalized by dividing by V_(s) (t=0) where t=0 is the instant at whichthe applied electric field is terminated. The normalized decay curve isthen fit with a set of N exponentials of the form: ##EQU3## where S_(i)is the fraction of total pore volume in which the time constant equalsT_(i). The number of exponentials fit to the IP decay curve depends onthe signal-to-noise of the IP decay curve. We have found that N between3 and 5 gives satisfactory results for typical signal-to-noiseconditions. FIG. 3 shows a typical example using four relaxation timesfor a Berea sandstone core. Each time constant is converted to a poresize using:

    A.sub.i =√DT.sub.i                                  (4)

where A_(i) is the radius of the i^(th) pore size. This results in Ndata pairs of the form (S_(i), A_(i)) from which a pore size histogramand cumulative pore volume versus pore size curve can be constructed.FIG. 4 shows the cumulative pore volume curve determined from the IPdecay curve derived from the aforementioned tool. If the additionalassumption is made that ##EQU4## where P_(c).sbsb.i is the capillarypressure in dynes/cm², γ is the interfacial tension between the wettingand non-wetting fluid (dynes/cm), θ is the contact angle, R_(i) is thepore throat radius in cm, A_(i) is the pore radius in cm, and h is aconstant relating pore throat radius to pore radius, then a capillarypressure curve can be constructed from the N data pairs (S_(i),P_(c).sbsb.i). This assumption is based on the relation between P_(c)and R for capillaries of circular cross section and by the assumptionthat pore throat radius and pore radius may be directly proportional toone another. The constant γ is typically 35 dynes/cm for oil/water, andcos θ=1 for water-wet rocks; γ=480 dynes/cm for mercury/air and cosθ=-0.7666, as is well known in the art. The constant h is typically inthe range 2.5-4 and may be determined more accurately from correlatingcapillary pressure measurements and IP decay curves on cores. FIG. 5shows a Hg/air and oil/water capillary pressure curve determined fromthe induced polarization voltage decay of FIG. 3 (using h=2.5), as wellas the actual Hg/air capillary pressure curve measured for this sample.

Several empirical relations have been shown to reflect how permeability,K, depends on pore size and pore distribution. The simplestapproximation is that:

    K≅<A.sub.i.sup.2 >                               (6)

where K is the permeability in millidarcies and <A_(i) ² > is the meansquare pore radius in square microns. A more accurate empirical ralationis:

    K≅cφ.sup.m* <A.sub.i.sup.2 >                 (7)

where c is a constant, φ is the porosity which may be determined from aporosity log, m* is the cementation factor which is typically in therange 1.5-2.5, and <A_(i) ² > is the mean square pore radius determinedfrom the pore size distribution. In this relation, if K is inmillidarcies, A_(i) in microns, and φ is in decimal, then the constantc=20. From the Waxman-Smits relation (M. H. Waxman, L. J. Smits,"Electrical Conductivity in Oil-Bearing Shaly Sands", Soc. PetroleumEngineers Journal, June 1968, pp 107-122): ##EQU5## where F* is theformation resistivity factor. The formation resistivity factor can bedetermined from the in-phase resistivity measured by the inducedpolarization logging tool if the formation salinity is known. Thus, ifC_(w) is the conductivity of the formation brine (mho-meters⁻¹), C_(I)is the in-phase conductivity (mho-meters⁻¹) measured by the IP loggingtool, then: ##EQU6##

As an example of the application of these equations, the permeabilitycomputed for the Berea sandstone of FIG. 3 is 477 millidarcies usingequation 7 and 604 millidarcies using equation 8 (φ=0.21, m*=1.78,F*=15.8). The measured brine permeability is 425 millidarcies at 1.0Molar NaCl.

Other permeability relations based on the shape of the capillarypressure curve are also well known in the art (W. R. Purcell, "CapillaryPressures--Their Measurement and the Calculation of PermeabilityTherefrom", Pet. Trans. AIME, February 1949, pp. 39-48; J. H. Thomeer,"Introduction of a Pore Geometrical Factor Defined by the CapillaryPressure Curve", Pet. Trans. AIME 219, 1960, pp 354--358). If core plugsare available, permeability measured on the core plugs can be used tocalibrate the pore scale determined from the IP log measurements ofrelaxation time. Once the correlation between core measurements andin-situ measurements is made, monograms or tables can be provided forobtaining permeability directly from the borehole measurements. Thiscorrelation may be programmed into a computer algorithm so that theresponse of the induced polarization logging tool is displayed directlyin permeability units.

When both water and oil are present in a water-set formation atcapillary equilibrium, the oil phase occupies the larger pores and thewater phase occupies the smaller pores. This results in a reduction inthe IP decay curves at the longer times, corresponding to the largerpores entered by the oil phase. FIG. 6 shows the more rapid decay atlong times for a Berea sandstone with 50% oil saturation relative to thesame core at 100% water saturation. If equations 6 or 7 are utilized foran oil-bearing formation, the reduced permeability calculated will beindicative of the relative permeability to water at the particular oilsaturation present in the formation.

From FIG. 6 it is evident that comparison of the IP decay curve for the100% water-saturated core and the IP decay curve for the oil-bearingcore can be used to determine the pore size entered by the oil phase.Thus, in FIG. 6 the oil phase has entered all pores larger than thepores corresponding to a relaxation time of 80 milliseconds. The 100%water-saturated core IP measurements can be obtained from core materialfrom the formation of interest which has been extracted of brine andhydrocarbons and resaturated with formation brine. In addition, ifcapillary pressure curves are measured on the cleaned core material, thein-situ oil saturation can be determined by comparing the 100%water-saturated core IP measurement with the logging IP measurements.The capillary pressures can be computed for the relaxation time usingequations (4) and (5). Entering the capillary pressure curve of FIG. 5at this capillary pressure gives the in-situ oil saturation and theheight of the oil column by well known techniques (E. J. Lynch,"Formation Evaluation", Harper and Row, New York, 1962, p. 43).

The method of the present invention also applies if a repetitive bipolartime-domain waveform or a bipolar time-domain waveform with dead time isutilized. As is well known in the art, a decay curve from these andsimilar waveforms can always be constructed by linear superposition ofthe basic step function response described above (J. R. Wait,"Overvoltage Research and Geophysical Applications", Pergamon Press, NewYork, 1959).

The use of a focused induced polarization logging tool has severaladvantages in accordance with the teachings of this invention. The IPdecay curve can be obtained with the electric current from theaforementioned IP logging tool focused either horizontally orvertically. According to the teachings of this invention, sinceelectrolyte diffusion occurs substantially in the direction of theapplied electric field, separate horizontal and vertical pore sizedistributions and permeabilities will be obtained. If a conventionalinduced polarization looping tool without focusing is utilized, separatevertical and horizontal values cannot be determined. In addition, thedepth of investigation, thin bed resolution, and borehole correctionsare all poorer with a non-focused IP logging tool. Nevertheless theother teachings of this invention can still be applied to logging with anon-focused IP logging apparatus.

An alternative method of the present invention is to operate the IPlogging tool in the frequency domain. This requires that the tool beoperated at several substantially different frequencies in the lowfrequency range, preferably between 10⁻³ and 10² Hz, and adjacent theformation of interest. The IP measurements at the different frequenciescan be made either sequentially or simultaneously by using asuperposition of the different discrete frequencies. The phase shift foreach frequency is measured and plotted to obtain a curve representingthe phase shift versus frequency over a substantial range offrequencies.

Referring to FIG. 7, there is shown a log-log plot of phase shift versusfrequency for three different formations. The formation A has thelargest pore dimensions, while the formation C has the smallest poredimensions. Point a in FIG. 7 shows the frequency at which the phaseshift in formation A starts to decrease at least as fast as f¹.0. Thisfrequency corresponds to the longest time constant in the IP decay curveand also to the largest pore size in the formation.

The phase shift versus frequency curve can be decomposed into a set ofamplitudes and time constants, using linear transform methods well knownin the art. The amplitudes and time constants are the same as those fromthe IP decay curve, and the method of this invention proceeds as aboveto determine pore size histogram, capillary pressure curve, andpermeability.

A shortened form of the above procedure is to obtain the phase shiftversus frequency plot as described above, and measure the frequency,f_(min), at which the phase shift starts to decrease substantially asf¹.0. This frequency f_(min) can be used in the following simpleapproximation for formation permeability: ##EQU7## where K is thepermeability in millidarcies, D is the diffusion coefficient in cm²/sec, and f_(min) is in Hz. Other empirical relations relatingpermeability to f_(min) are also possible as derived from correlationsbetween measured permeability on core plugs and f_(min) determined fromthe IP log. This approach is illustrated in FIG. 7 wherein the phaseversus frequency curves for three formations A, B and C are shown. Thef_(min) for the three formations are indicated at a, b and c,respectively.

What is claimed is:
 1. A method for the in-situ determination ofpermeability of a formation, comprising the steps of:logging theformation of interest with a time dependent induced polarization loggingtool to produce transient induced polarization measurements; convertingthe transient induced polarization measurements obtained by the loggingtool to an induced polarization formation response function; determiningthe pore size distribution from the induced polarization responsefunction; and determining the permeability from the pore sizedistribution.
 2. The method of claim 1 and in addition determining thecapillary pressure histogram from the pore size distribution.
 3. Themethod of claim 1 wherein said logging step consists of inducedpolarization measurements made in the time domain, and said responsefunction is a time domain voltage response function.
 4. The method ofclaim 3 wherein said logging step consists of induced polarizationmeasurements made using a current-on time and current-off time in therange from one 1/100 to 100 seconds duration.
 5. The method of claim 4wherein said induced polarization (IP) time domain voltage responsefunction is normalized by dividing by the IP decay voltage at the timethe current flow is terminated, and the normalized response function isthen fit with N exponentials in the expression: ##EQU8## where V_(s) (t)is the IP decay voltage at time t, V_(s) (t=0) is the IP decay voltageat the time after the primary charging voltage is terminated and S_(i)is the fraction of pore volume whose time constant equals T_(i).
 6. Themethod of claim 5 wherein a pore size distribution is determined fromthe equation:

    A.sub.i =√DT.sub.i

where A_(i) is the radius of the i^(th) pore size of pore volumefraction S_(i) for which the time constant equals T_(i), and D is thediffusion constant of electrolyte ions in the formation brine in cm²/sec at formation temperature.
 7. The method of claim 6 wherein acapillary pressure histogram consisting of data pairs (P_(c).sbsb.i,S_(i)) is derived from the equations: ##EQU9## where P_(c).sbsb.i is thecapillary pressure in dynes/cm² for the i^(th) pore throat size, S_(i)is the pore volume fraction corresponding to P_(c).sbsb.i, γ is theinterfacial tension between wetting and non-wetting fluids in dynes/cmat formation conditions, θ is the contact angle, R_(i) is the porethroat radius in cm, and h is a constant relating pore throat radius topore radius.
 8. The method of claim 6 wherein the permeability K isderived from the equation:

    K≅<A.sub.i.sup.2 >

where K is the permeability, and <A_(i) ² > is the mean square poreradius.
 9. The method of claim 6 wherein the permeability of K isderived from the equation:

    K≅cφ.sup.m* <A.sub.i.sup.2 >

where K is the permeability, c is a constant, φ is the fractionalporosity, m* is the cementation exponent having a value between 1.5 and2.5, and <A_(i) ² > is the mean square pore radius.
 10. The method ofclaim 9 wherein logging measurements made with resistivity and densitylogging tools opposite the aforesaid formation are used to obtain thefactor φ^(m*).
 11. The method of claim 1 wherein said logging stepconsists of induced polarization measurements made in the frequencydomain, at a plurality of frequencies, and said response function is arecord of induced polarization phase shift versus frequency.
 12. Themethod of claim 11 wherein said logging step consists of inducedpolarization frequency domain measurements made in the frequency range0.001 Hz to 100 Hz.
 13. The method of claim 11 wherein the permeabilityis determined from the following equation: ##EQU10## where K is thepermeability, D is the diffusion constant for electrolyte ions in theformation brine at formation temperature, and f_(min) is the frequencyat which the frequency response function decreases substantially asf¹.0.
 14. The method of claim 11 wherein the frequency response functionis decomposed into a set of amplitudes and time constants using lineartransforms and the resulting amplitudes and time constants are used inthe method of claim
 5. 15. The method of claim 11 wherein the inducedpolarization frequency domain measurements are made simultaneously at aplurality of frequencies using a waveform which is a superposition ofthe desired frequencies.
 16. A method for determining in-situ thepermeability of a formation surrounding a borehole, comprising:loggingthe borehole with an induced polarization logging tool having at leastone focused source electrode; during said logging step applying currentfrom said at least one focused source electrode to the borehole at aplurality of frequencies in the range 0.001 Hz to 100 Hz while thelogging tool is opposite the formation of interest; measuring the phaseshift of the electric field produced by said current at the variousfrequencies; transforming the phase shift versus frequency measurementsinto amplitude versus time; determining the pore size distribution for Nsizes from the expression:

    A.sub.i =√DT.sub.i, i=1 to N

wherein A_(i) is the radius of the i^(th) pore size, T_(i) is the timeconstant for the i^(th) pore size, and D is the electrolyte diffusionconstant at formation temperature; determining a capillary pressurehistogram consisting of data pairs (P_(c).sbsb.i, S_(i)) from the poresize distribution using the expression: ##EQU11## wherein P_(c).sbsb.iis the capillary pressure, γ is the interfacial tension between wettingand non-wetting fluids, θ is the contact angle, R_(i) is the pore throatradius, and h is a constant relating pore throat radius to pore radius;and determining the permeability from the pore size distribution usingthe expression:

    K=cφ.sup.m* <A.sub.i.sup.2 >

where c is a constant, φ is the fractional porosity, m* is thecementation factor, and <A_(i) ² > is the mean square pore size.
 17. Amethod for determining in-situ the permeability of a formationsurrounding a borehole, comprising:logging the borehole with an inducedpolarization logging tool having at least one focused source electrodewhich emits current into said formation; measuring during said loggingstep the phase shift of the electric field produced by said current at aplurality of frequencies in the range 0.001 Hz to 100 Hz; determiningthe highest frequency f_(min) at which the phase shift decreasessubstantially as f¹.0 ; and determining the permeability from thefrequency f_(min).
 18. The method of claim 17 wherein the permeabilityis determined from the expression: ##EQU12## where K is thepermeability, D is the diffusion constant for electrolyte ions in theformation brine at formation temperature, and f_(min) is the frequencyat which the frequency response function decreases substantially asf¹.0.
 19. The method of claim 17 wherein the focused electrode array isadjusted to obtain separate vertical and horizontal permeabilities. 20.The method of claim 17 wherein the permeability is determined using anomogram relating permeability and f_(min).
 21. The method of claim 1wherein said logging step consists of measurements made using a focusedelectrode induced polarization logging tool, containing at least onefocused source electrode.
 22. The method of claim 21 wherein saidfocused electrode logging tool is focused alternately in a horizontaland vertical direction, and the permeability is determined separately inthe horizontal and vertical directions.
 23. The method of claim 1wherein the induced polarization respons4 function is correlated withvalues of permeability obtained from core measurements to produce adirect output permeability log.
 24. The method of claim 1 wherein theinduced polarization response function is correlated with capillarypressure curves obtained from core measurements to produce a directoutput capillary pressure curve.
 25. The method of claim 17 wherein theinduced polarization response function is correlated with values ofpermeability obtained from core measurements to produce a directionoutput permeability log.
 26. The method of claim 17 wherein the inducedpolarization response function is correlated with capillary pressurecurves obtained from core measurements to produce a direct outputcapillary pressure curve.